Ensembles of physical states and random quantum circuits on graphs 
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In this paper we continue and extend the investigations of the ensembles of random physical states introduced 
in A. Hamma et al. |Phys. Rev. Lett. 109, 040502 (2012) |. These ensembles are constructed by finite- 
length random quantum circuits (RQC) acting on the (hyper)edges of an underlying (hyper)graph structure. The 
latter encodes for the locality structure associated with finite-time quantum evolutions generated by physical 
i.e., local, Hamiltonians. Our goal is to analyze physical properties of typical states in these ensembles, in 
particular here we focus on proxies of quantum entanglement as purity and a-Renyi entropies. The problem 
is formulated in terms of matrix elements of superoperators which depend on the graph structure, choice of 
probability measure over the local unitaries and circuit length. In the a = 2 case these superoperators act on 
a restricted multi-qubit space generated by permutation operators associated to the subsets of vertices of the 
£N| graph. For permutationally invariant interactions the dynamics can be further restricted to an exponentially 

smaller subspace. We consider different families of RQCs and study their typical entanglement properties for 
finite-time as well as their asymptotic behavior. We find that area law holds in average and that the volume law 
is a typical property (that is, it holds in average and the fluctuations around the average are vanishing for the 
large system) of physical states. The area law arises when the evolution time is O(l) with respect to the size L 
\Q of the system, while the volume law arises as typical when the evolution time scales like O(L). 
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I. INTRODUCTION 

The study of the statistical properties of ensembles of pure quantum states is an important topic in quantum information theory, 
quantum statistical mechanics, and quantum many-body theory. The ensemble of pure quantum states can be chosen to mimic 
the uniform distribution of states in the Hilbert space, of low energy states of random Hamiltonians, or states that can be obtained 
by some random quantum evolution. One can combine group theoretical and statistical tools by construction the ensemble with 
group theoretic-methods. An important example in quantum information theory is the use of ensembles of random unitary 
operators to perform quantum algorithms. In this case one picks the unitaries from the Haar measure on the unitary group. 
A related ensemble is the ensemble of states in the Hilbert space that can be obtained by some random preparation. If one 
is allowed to obtain all the possible states with the same probability, one has again used the Haar measure over the ensemble 
of states. Recently, this kind of ensembles has been studied in relation to questions of typicality of the expectation value of 
observables and the foundations of statistical mechanics 1 1 1-[ 14| 

In this paper, we engineer and analyze ensembles £ of pure quantum states for multi-partite systems that incorporate two key 
features in quantum information theory: randomness and locality. In practice, £ is generated by applying a Random Quantum 
Circuit (RQC) to a reference state. The RQC is obtained by picking stochastically subspaces of the total Hilbert space that obey 
some locality constraint given by an underlying graph-theoretic structure. These subspaces serve then as the support of random 
unitaries. In this model, randomness enters two times, in the choice of the support of the unitary and in the choice of the unitary 
itself. The ensemble so generated finds physical motivation in approximating the evolution of a multi-partite system given by a 
(time dependent) local Hamiltonian |[T5l[T6l . Our analysis is conducted by putting together graph theoretic, group theoretic and 
operator algebra tools. In particular, we show how to encode all the relevant information about the RQC in the action of a single 
superoperator. 

From many perspectives, an extremely important property of pure states of multi-partite quantum systems is their entangle- 
ment fT7l[T8l (in a given bipartition AUBof the system of linear sizes La and Lb)- To measure entanglement in the ensemble 
£ we study the a— Renyi entropies of the reduced density matrix to the subsystem A. In particular, because the linear entropy 
is a lower bound for the Von Neumann entropy, the case a — 2 is very interesting. We show in the following that the average 
purity in the RQC defined ensembles attains asymptotically the minimum possible value, and therefore gives a tight bound to 
the Von Neumann entropy. 

As we shall see in the following, we are interested in studying two different regimes. The first regime is obtained when the 
RQC is applied in one shot, or in a number O(l) of instances. The second regime, is obtained by applying the RQC a number 
of time scaling with the size La of the subsystem. 

We analyze two models: (i) the nearest neighbour Random Edge Model (REM) picks the support of the circuit in subsystems 
corresponding to the edges of a graph. This model leads to estimates about the time evolution of Renyi entropy when the system 
evolves with Ultra-local Hamiltonians. A variant of the REM, which we call the Fully Connected REM, where the local structure 
is given by a completely connected graph, elucidates the utility of the superoperator approach for purity dynamics by mapping 
the calculations to an exponentially smaller space, (ii) the second model of interest is the Contiguous Edge Model (CEM). In 
this model, we again specify a graph and then we consider a random quantum circuit which has support on the all N qubits of 
the system (and therefore the depth of the circuit is N). The random circuit is thought factorizable into k terms that only have 
support on the edges of the graph. The application of such random circuit to a reference state will here be referred to as a cycle. 
Cycles can be iterated a number n c of times. The CEM is intended to produce ensembles of states that come from the unitary 
evolution induced by a local Hamiltonian starting from completely disentangled fiducial states <f>. In this case, n c represents 
discretized time. For regular graphs of linear size L, we look at the reduced state pa and we show in the CEM typicality of Area 
law of Entanglement for small times, n c = 0(1) and Volume law for times n c = O(La)- The calculations of the Purity lead 
to an instructive algebra of the swap operators on different subsets A of the total nodes in the graph T = (V, E). Moreover, 
asymptotically in n c we show that the reduced system becomes the completely mixed state. 

Note that in the cases where the RQC is applied a number O(La) of times, our result shows that the ensemble attains in 
average a reduced system which is close to the completely mixed state. This property is shared with the ensemble of states 
over the Haar measure. In this sense, we show that the ensemble £, even if contains a small fraction of the states in the Hilbert 
space, nevertheless, locally realizes the averages over the Haar measure. This concept is familiar in the context of t— designs. 
A t— design is an ensemble {pi,ipi} of states that reproduces moments over the Haar measure, that is, Y^iPidi*) (V 7 !)®* = 
L(|V>) (tpD^dip. A unitary t- design is similarly an ensemble of unitaries such that £\ piU^piUj)^ 1 = j v U® 1 p(W) 9t dU '. 
In other terms, the average with the Haar measure of any polynomial function of degree t can be obtained by a t— design. 

It is known that exact t— designs require an exponential number of states (or unitaries). On the other hand, approximate 
t— designs are much more easily obtained. Indeed, it has been shown in the literature [ 19 1-[32 1 that RQC of size n 2 (where n is 
the number of qubits in the system) are approximate t— designs, for t = 1, 2, 3. The typical scheme for such RQC is to consider 
a random circuit where two qubit-unitaries are drawn with the Haar measure on U(A) for every pair of qubit in the system. 
Recently, it has been shown that even with some locality constraint RQC of linear size in n are up to 3— designs [33). 

It is important to emphasize that in this work, we focus on a particular problem, which allows us to obtain results beyond the 
asymptotic case. We are not trying to approximate any polynomial that is function of ip over the Haar measure, but, given a 
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bipartition AU B of the system, only the reduced system pa of linear size La- With this restriction, we find very powerful tools. 
In particular, (i) we are able to make statements about typicality of entanglement for circuits of any depth. Even the asymptotic 
case, scales with La being the size of the subsystem, (ii) Our results are valid for every momentum t of the statistical distribution. 
While the REM is a RQC with locality constraints, the protocol of the CEM is quite different, as we will see. It is motivated and 
inspired by the evolution induced by a distribution of time-dependent Hamiltonians. For this reason, the ensembles produced by 
the CEM considered here are of physical relevance for applications in the foundations of statistical mechanics. 

On the other hand, when the RQC is applied a O(l) number of times, we obtain an ensemble of states with typical area law 
for the entanglement. In some sense, this ensemble shares a lot with the set of the ground states of local Hamiltonians (without 
topological order). Indeed, all the ground states of such Hamiltonians can be obtained by a quantum circuit of fixed depth from 
some completely factorized state ll34ll . Whether or not in two spatial dimensions there is an area law for the entanglement in 
gapped systems is an important open problem in quantum many body physics [35 1. Our approach shows, that in ensembles that 
contain such ground states, the area law is typical. It is known that such states have the area law as upper bound. Indeed, the 
technique of the Lieb-Robinson bounds has shown that entanglement that can be produced in a subsystem A by evolving for a 
time t with a local Hamiltonian is upper bounded by 0(|cL4|i) [36], [37]. Our study shows that such upper bounds are saturated 
in average, and that the fluctuations are small. 

The structure of the paper is as follows: In section (ii) we describe the setup for our models, In section (iii) we briefly review 
essential mathematics for the remainder of the paper, In section(iv) and (v) we present detailed studies of two particular models 
and conclude in section (vi). 



II. SETUP 

Our scheme to investigate typicality of entanglement involves two elements: 

1. A (hyper)graph T = (V,E) whose nodes V represent local Hilbert spaces corresponding to local degrees of freedom of 
the multi-partite system and where the edges E represent the support of interactions. 

2. A Random Quantum Circuit (RQC) that acts on an input fiducial state <f> that can be represented conveniently using the 
abstract tensor product Hilbert Space on the nodes of the Graph. 

The system S is defined by the tensor product of many local Hilbert spaces H — <E) V £vH v . We are concerned with finite 
d— dimensional Hilbert spaces H v , or qudits. The total number of qudits in the sistem is N — \V\. We can regard the set V 
as the set of vertices of a (hyper)graph T = (V,E). We remind that a hypergraph is a set of vertices with a collection E of 
subsets of V called edges. Mathematically, E is a subset of the power set "P(V)\0. For instance, the usual graph, is a set with 
a collection of pairs. The hyper graph is a natural structure for multi-partite quantum systems because we can associate the 
vertices to the Hilbert spaces of each particle and the edges to the support of interaction terms in the (for instance) Hamiltonian. 
A bipartition in S is introduced by considering a bipartition of the set of vertices V — A\J B and then considering the tensor 
product H — Ha ® Hb, where Hy = ®x£Y<zvH x . If c? = dim "Hz for every x, then dim'Hy = e?' r '. If X g E is an edge of 
the (hyper)graph, X = {xi, x\x\} with Xi £ V and |X| is the cardinality of the edge X. We define the Hilbert space with 
support on X as Hx — H Xl ® . . . <£> H x , x , ■ 

One can regard the subsystem A as the physical "system" of interest, and its complement B as its "environment", usually 
assuming that dA <C ds- A totally factorized state for the whole S (system +environment) can thus be written as \<p) — 
®i£A \<t>i) ®i£B \4>j) where obviously G Hi- In this work, the fiducial state is any totally factorized state (and which 
exactly does not matter, as we shall see). 

The fiducial state is the input to the RQC that picks edges X in E to act on, according to some probability distribution 
V{X), with a unitary operator U acting on Hx- The unitary Ux is picked with some measure dp(U\X), e.g. the Haar measure 
over U(d) x]l ). In other words, we first pick an edge X g E with probability V(X) and then we pick a random (with measure 
dp) unitary with support on Hx ■ The RQC itself can thus be labelled by the probability distribution of the edges and choice of 
the measure over unitaries. The ensemble £ is then completely specified by varying on the fiducial states <\> and the chosen RQC, 
in the following way 

£(V,dp) = {U x \<f>)}x,u x ,<j> (1) 

More generally, one can describe a general Random Quantum Circuit where instead of picking just a single edge at each 
step the circuit chooses a subset of the nodes. Such RQC can be described using the joint probability distribution : 
(Xk,Xk-i,...,Xi) -> V&^X^Xk-u-.tXi) € [0,1] where I, C y,i = 1,2, ...,k are subsets of the set of vertices V 
of the (hyper)graph, or, in other words, edges in E. The sequence of set valued Random Variables, X\, X2, Xk, can 
be seen as a stochastic process of length l k' . Assuming that such a selection of subsets of V is a Markovian process, one 
can express the action of RQCs of arbitrary depth using just the Markov Transition matrix, (Xk\Xk-i) which satisfy 
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J2 Xa MW(X a \Xp) = IMi = 1,2,..., k where a, = 0, l,...,2l y l - 1 label the elements of the powerset of V. 
The joint probability of choosing the set X, at the Yth step where T ranges from i = 1,2, k is given by 

v^{x k ,x k _ 1 ,....,x 2 ,x 1 ) 

= M^(X k \X k _ 1 )V^ k - 1 \X k _ 1 ,X k _ 2 ,....,X 2 ,X 1 ) (2) 

where V^ k ~^ (X k _Xi - ,^1) = J2x k i^k, X k -i, is the marginal distribution of the preceding k — 1 set-valued 

random variables. Iterating this equation one can write the joint probability at the fc'th step as a product of the transition matrices 
given an initial probability vector (Xi) as follows: 

7>(*>(X fc) X fc _ 1 ,...,Xi) 

= M( fe )(X fe |X fe _ 1 )M( fc - 1 )(X fe _ 1 |X fe _ 2 ) x 

x M {2 \X 2 \X 1 )V {1) {X l ) (3) 

where (Xi) is a column vector representing the probability of choosing X\ C V in the first draw by the RQC. In particular, 
for independent choices of Xi , i = 1, .., k at each level of the circuit, we have that the elements of the Markov transition matrix 
satisfy, AfW(Jr fc |Jr fc _i) = V®(X k ), which implies V (k) (X k , X k _ x , ....,X X ) = ]J k =1 V {i) {Xi). The depth of this circuit is 

Yli—i If, as we assume, the subsets Xi are finite, the depth of this circuit is O(k). 
The associated fc— iterated ensemble is then given by 

= {U Xk ....U X2 U Xl \$)}x,u x ,$ (4) 

where dfXi, i = 1, k are the measures with which the RQC chooses the unitaries to act on the corresponding chosen subsets 
Xi C V. 



III. ENSEMBLE STATISTICAL MOMENTS AND SUPEROPERATOR FORMULATION 

Once the ensemble £ k of physical states has been constructed, we need to derive the associated ensemble for the reduced 
system A. This is naturally obtained by tracing out the environment B, that is: 

£ k A = {T lBP \ P e£ k }. (5) 

In order to compute the Renyi entropies for this ensemble, given a density matrix 6 £\, we will compute the trace of its 
a— power 

p a (n) = TY A (n a ) (6) 

with a > 1. The central objects of our analysis are the statistical moments of the P a (Q.) within £\, for instance the average 

P a (fl) and higher moments. We want to stress that the knowledge of all the statistical moments in £\, is equivalent to being 
able to compute also all the statistical moments of observables with support on the reduced system A or polynomial functions 
/(f2) of arbitrary degree. One pertinent example is that knowledge of the a—powers allow us to compute the Renyi entropies 
H a := j3^1og[Tr(p^)], which are continuous w.r.t the parameter a [38 1. Also, < and therefore H a= \ > H a=2 . Since 

the VonNeumann entropy E(p A ) can be obtained as the limit, E(p A ) ■— \ira a ^\+H a we see that the 2-Renyi entropy lower 
bounds the VonNeumann entropy. Moreover, for close to minimal purity the bound gets very tight 1 39 1 . In order to compute the 
statistical moments in £\, we will use quantum information theoretic tools, then we will introduce a superoperator formulation 
which will allow for a compact description of the statistical properties of the reduced system, and will show how locality of the 
interactions influences its entanglement properties. We start by recalling that for every density matrix f2 of the reduced system 
A, 

P a = TY A (fl a ) = Tr A (n® a f A a) ) 

= Tr[uj® a T A a) ] (7) 

where : (H A )® a n- (%^)® a , \i\,i 2 , —.,i a ) \i a , H, ■ — , i a -i) is the order a shift operator acting on the Ha 

subspace alone. It is thus arestriction of the operator T A a ^ to just the A subspace where \tn A )®a = and T A {(-Hb) 9 " = 
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1(^ B )®<*. Here oj is the state of the total space f2 = Ti'b[cj]. Note that up to a rearrangement of spaces = <S)ieAT} a <8> Is. 
To avoid later confusion we remark at this point that we use tensored copies '%® Q of the total Hilbert space H = ® v ^vT~Lv that is 
itself a tensor product of local spaces. We denote the former as rank-a tensor space while the latter just as a tensor product space. 
Note that a state in H® a of the form ift® a is symmetric under permutation operators acting on the a , H spaces. In particular 
the a-tensored copy of the totally factorized state uj® a is symmetric under T^* , A <Z V. Starting with a fiducial (completely 
factorized) state lj = \(f>) (<fi\ = <E>i£v (</>i|> tn e average of P a over unitaries U x that act on the set of vertices X drawn by 
the RQC is given by 

W Ux = J dn(U\X)Tr[(U x ujU x )® a T { A a) ] 

= J dii{U\X)Tr[uj® a {U x )® a T^U® 01 } (8) 

where in going to the second line above we used cyclicity of the trace. By choosing the Haar measure as d[i in the integral above, 
we can perform integration using standard results from Representation theory of groups Il40l - ll42ll . 



A. General formulation 



In this section, we present a formulation in terms of superoperators in order to treat general RQCs. The key insight is given 
by examining the action of the operator valued integral, TZ X '■ B{%y a ) — > B(Hy a ), defined by 

nP(d) = [ d^{u x ){u x r a 6uT (9) 



where as before U x are unitaries drawn with a measure d/i(U\X) acting on the set of vertices X drawn by the RQC and 
O € B{'Hy a ) is any operator on the same space. We are interested in the specific case where O = Ti . Note that TZ^ * s 
not necessarily a self-dual (hermitian) superoperator and whether it is so depends on the choice of unitaries and the measure of 
integration. The stochastic procedure to pick the edge X is encoded by defining the superoperator 

n {a \d) = E v{x)n { x \6). (io) 

XcV 



Equation ( 10 1 represents the action of a RQC of depth O(l). For a circuit of arbitrary depth O(k), the superoperator takes the 
form: 

A;. ACV 

i — k 

Y[dn(U\Xi)(U Xk U Xl )^ a d(U Xk U Xl ) 9a (11) 



In the case of uncorrected choices of the subsets at each step eq.( 1 1 1 takes the form: 

<'(o)= e n^w* 

X,. ,X k CVi=l 

i — k 

Y[df,(u\x t )(u Xk u Xl )^ a d(u Xk u Xl T a (12) 

From eq.(JT^j» one can see that for the same case (of independent choices of Xi), (O) — Yii=i T^i** O n the other hand 
if the choices of X t are fully correlated i.e. V {i) = V , V i = 1, .., k then Tl^UO) = {K^) k (6). 

The superoperator formulation of the averaging procedure over the unitaries allows one to write the average a-moment for a k 
level RQC very compactly as (c.f. Eq(8| : 

p^H^^SW)) (B) 

where (Oi, O2) = Tr(Oj02) is the Hilbert-Schmidt inner product. 
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B. Single Edge Model 

As an illustrative example, we show how the averages Eq.([8]l can be obtained for the simplest model consisting of two 
subsystems Ha = ® r H r and Hb = ® S H S connected by a single edge. The graph consists of partitions A and B connected 
by an edge with V = A U B and E = {{i, j}}, and i G A, and j £ B. The 2-body unitary Ux has support on the edge 
X = {i, j}, that is, the Hilbert space Hx = Hi ® Hj, where dim 7^ = dimHj = d. We call this single edge X. For sake of 
simplicity, we show first the calculation for a = 2 and generalize to any a in the following. We have 

'X J 



pa =2 u =T r[H® 2 J dU x {U x f 2 T A 2) Uj 

= Tr[H^ 2 J dUx(U x fX i2) U r fT A %] 



Trl^r* N d (T A %+T A %)} 



= N d Tr[(u/)® 2 T^.} +7V d Tr[H^ 2 T^.] = 27V d (14) 

where Nd — (,p+i) an( ^ ^\ is the dimension of the Irreducible subspace labelled by A = ±, and <S)\H\ = (Hv)® a - The 
II + , II_ are the projectors onto the totally symmetric and totally anti-symmetric subspaces of (Hi ® Hj)® 2 . In the third line of 

(2) 1. -T (2) T< 2) 

the above derivation the trace of the product of the operator Tj> ; with the projector II_, is zero since II = — - — ^ — 1 — => 

(2) ('2\ (2) 

Tr(T^ II_) = (l/2)TTi j[T^ 'tj — IjTL- ] = using the fact that the swap operators on the same subspace square to one i.e. 



(Tjj) 2 = tij. In order to get to the last line in eq.( 14 1 above we use the fact that both the traces yield 1 since the T x ' , X (Z V 
operators acting on the totally symmetric state (cj)® T (under permutations of the rank-2 tensor space H® 2 ) leave it invariant. 

Note that the fact that Ux has support on both Ha andHs is crucial. Indeed, if we consider a unitary U a with support on 
just A (and similarly on E), we would obtain / dUTv[u® 2 (U A )® 2 T A 2) Uf 2 ] = 1 since T {2) E S 2 (symmetric group of order 2) 

is also an element of the commutant of U® 2 i.e. for U — Ua,b, [U a 2 b , T a ^] = 0. This property will constitute the basis for the 
general formulation to describe general RQC. 

This calculation can be generalized to all a > 3. The average of P a over Ux € U(Hx = Hi <£> Hj) is: 

W = Tr[io® a / dU(U<<)® a T A a) U® a } 



(2) 



= Tr[^-^±Tr[TllT A a) ]] = Tr[^-^±]TT[Il%T A a) } 

a_[_ 

a -\ d 2 -l ) aeS a 
I d?-\ ) a£S a 

,2 



1 x - m x 

V d 2 -l / A:IRR of S a 



(15) 



where II" and (d+) are the projector and the dimension of the totally symmetric part of Hf a , respectively. m\, d\ are the 
multiplicities of the different Irreducible representations (IRRs) of S a in the 'a' representation. Xa(^i) is the character in 
those IRRs of cr, which is the fixed operator that performs the action Hf a — > Hf a by taking a basis state i 2l i a ) — > 

\i<xi ill H) i •••) ia-i in other words the T A a ^ operator restricted to just the i'th subspace. In going to the second line above we 
use the fact that the trace of u>® a w.r.t. all other projectors except onto the totally symmetric one is zero. To prove the last line 
we use Schur's Orthogonality theorem flOl (or see appendix [A). 
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Being able to compute averages of arbitrary powers of the reduced system, we can also compute higher statistical moments. 
We are especially interested in the variance of the purity: 

Var{P a = 2 ) = (P«= 2 ) 2 - (P^ 2 ) 2 (16) 

For the Single Edge Model the first term is computed as follows 

(P«= 2 ) 2 = J [dC/ e ]Tr[^ 4 Vf 4 T { 2 2) <g> T^ 34) Uf 4 ] 

= Tr[^±]Tr[n + (T l (12) ®Tf 4) )] (17) 
«+ 

where d + = C^+f _1 ) = is the dimension of the totally symmetric subspace of (Hi <8> 'Hj)® 4 and 

II + = J, (X) CT eS4 °« ® a 'j) m trie projector onto it. In the above, only the symmetric projector is relevant because \<f> n ) (</>„|® 4 
has support only on it. The operators T^ 12 ^ , T^ 34 ^ represent swaps on the first two and the last two spaces again of (Hi ® Hj)® A 
i.e. in cycle notation the tensor product J^ 12 -* ® y\( 34 ) represents the (12) (34) element of S4. 
Direct calculations (see appendix 



B|) leads to (P«= 2 ) 2 = ^^(d^VA-i) 1 Finall y we obtain 



^ [pQ " ] - ( rf 2 + 3)( ( ' 2 +2 T( rf 2 + i) 2 (18) 



If the subsystems i, j are qubits, d — 2 and we obtain Far [P" ] = .017 



We have seen from eq.( 14 1 that for a = 2 the superoperator on a two-vertex edge X := {a, 6} and any subset A C V of the 
graph acts like 

nf{T A ) = N d (T AXX + t Aux ) ini^AinB^ 

n x r) (T A )=T A otherwise (19) 

From this key relation we see that the X— supported superoperator takes a permutation operator T A with support in A and yields 
two permutation operators with support in A\X and A U X, respectively, when X is across the boundary between A and B, 
otherwise is trivial. This action has two implications: (i) a reduction of purity (because the number Nd < 1 appears). Therefore 
only when IZx is supported across the boundary it can have a non trivial effect, and (ii) a propagation in the bulks of both A and 
B at a distance given by the diameter of X. When iterating with k, this implies a propagation of entanglement in the bulk that 
scales with k, as we will see in the concrete models described in the next sections. In the case of general hyper graphs, the edges 
may contain any number of vertices. 

In the following, we will exploit the power of this formulation in order to study models defined on non trivial graphs. In this 
way, we will be able to show typicality of area law and volume law for the entanglement of a subsystem. 



IV. RANDOM EDGE MODEL 



A. Random Edge Model for a general graph 

In this model, we consider a graph (V, E) where the set of edges E C V 2 , that is, a usual graph. The RQC draws an edge 
X G E according to the flat unit normalized measure P:Ic£-) V{X) € [0, 1], with V{X) = r|r. Conditioned to the 
extraction of the edge X, a unitary with support on X is drawn with the Haar measure d[iHaar(U\X). In this way, we obtain 
the ensemble 

£CP = ,dflHaar) = {U X \^)} X ,U X ,<P ( 20 ) 

where as usual (f> is the completely factorized fiducial state. This ensemble is obtained by extracting RQC of depth < 
m&xxeE \X\ — O(l). If we denote by dA C E those edges that go across the boundary, that, is those that have non-null 
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intersection with both A and B, the average purity is given by 



pa=2 = V(X)Tr[cu® 2 j d[U x }(U x )® 2 T A (U x f 2 } 



XeE\dA 

+ Y, V(X)Tr[u® 2 j dU x {U x r 2 T A {U x r 2 ] 

XedA J 

= y nx)xi 

XeE\dA 

+ ]T V(X)[N d TY[^ 2 T AuX ]+N d TT[oj^ 2 T A \ X }} 

XedA 

= (l-q) + 2qN d (21) 

where q = ^2 X edA = j^j is the net probability of a boundary node of A interacting with one in B. Since local unitaries 

completely internal to A or B i.e. with no support across the bipartition cannot affect the purity we get a 1 for the integral in 
the first term on the R.H.S.. Note that this averaging procedure over the edge unitary executed once leads to an algebra of the 
Tx operators where with probability q — jA me RQC generates an equal superposition of T Aue and T A \ e (e being the edge on 
which the unitary acts) and with the complement of the probability (1 — q) it leaves the operator T A invariant. At this point, we 
are ready to show how to compute the variance of the purity P a = 2 in the above ensemble. 



By using the technique leading to Eq.( 18 1 and the probability distribution V(X) of the REM, We obtain (with / = (P Q - 2 )2 



defined in eq.( 17 1): 

Var[P^ 2 M ] = (l-q + q{P»= 2 ) 2 ) -{l-q + q2N d ) 2 
nq(l + I-4N d ), (for<z«l) 
_\dA\ 4d 5 -4d 4 +20d 3 - 18d 2 +24d-2 

~W (d 2 + 3)(d 2 + 2)(d 2 + l) } ( } 

This is the variance obtained by the distribution obtained for a single edge extraction. We want now to compute the statistical 
momenta for the ensemble in which we extract the edges a number k of times. We iterate the procedure by using sequential 
independent identical RQCs and build the k— iterated ensemble £( k \ Let us consider unitaries of the form U — Yli=i Ux t 
where the JQ's and t/j^'s are drawn according to the probability distributions ]X=i ^C^*) anc ' IliLi d^iUxi respectively. If 
we now assume that the degree of each vertex is o(|<9A|) which implies that the boundary length changes negligibly due to the 



algebra of eq.(21 1 i.e we may take \dA U e\ ~ \dA\e\ ~ |cL4| then for the second iteration of the procedure we find that the 



average purity is given by : 

P^= Y nX)[V{X')T A + V{X')N d {T AUe ,+T AV ,)] 

X£E\8A 
X'eE\dA 

+ Y V(X)N d [V(X')T AUe +V(X')T AV ] 

XedA 
X'eE\d(AUe) 

+ Y V{X')N d [N d {T AueUe , +T AUe \ e ,) 

Xed(AUe) 
X'ed(A\e) 

+ V(X')N d (T AVue , +T AVV ,)} 
~ ((1 -q) + 2qN d ) 2 (23) 

It can be shown that under the same assumption (of negligible change of boundary length) the purity for k iterations goes as 
P£ =2 = (1 — q(l — 2N d )) k . One can understand the physical content of these calculations by considering the thermodynamic 
limit of large \E\ >> \dA\ => q small. In this limit the average 2-Renyi Entropy H 2 can be lower bounded (using concavity 
of the log function) by the logarithm (base 2) of the average purity i.e.: 



H 2 := -logP fc Q=2 > -logP^= 2 

H 2 > -Mog(l - g(l - 2N d )) ~ qk(l - 2N d ) (24) 
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The number k of iterations corresponds to the time in our scheme (and in the circuit model), so Eq.([24| implies linear increase 
of the entropy with time. It also implies that the entropy is proportional to the boundary of the region \dA\. For this model, 
the a— Renyi entropy turns out to have the form H a > y^log[Tr[.P Q ]] oc \dA\ which can be seen by considering P( Q ) = 

(I - q) + q Tr[u® a f d[U x ](U x )® 2 T A a) (Ux)® 2 } = 1 + q(C(a, d) - 1) where C{a, d) = ^gj^l is a group theoretic factor 

less than 1 with IT + the projector onto the totally symmetric subspace of H® a . 

All of the above calculations can be reformulated in terms of the expectation value of the stochastic Hermitian superoperator 
K = J2xev VWR-x = J2xev V i X ) I dU x (U r )® 2 T x U® 2 . In the next sections, we show that the asymptotic analysis of 
1Z for this case leads to a asymptotic purity of the form as in Eq.([24)>. 



B. Random Edge Model for the complete graph A'jv 

Let us now consider the case in which the graph of the REM is the complete graph Km, that is, a graph in which every two 
vertices are connected by an edge: E = V x V. The number of edges is of course 

\ E \ = \V 2 \ = (2)- As before, the RQC 

draws an edge X according to flat unit normalized measure V : X e E -)• V{X) e [0, 1], with V{X) = ^ . Conditioned to the 
extraction of the edge X, a unitary with support on X is drawn with the Haar measure dfj,jjaar(U\X). The analysis of this case 
highlights the utility of the superoperator formalism. We introduce a bipartition in the system by V = A U B with A n B = 0, 
with cardinalities \ A\ = Na and \B\ — Nb, and of course N = Na + Nb- The average purity after the k-th iteration is given 
by: 

^2(T A ) = Tr[w® 2 74 2) (Tf 2 )] (25) 

To see what form 7Z takes we have to see that after each extraction, if the extracted edge is not straddling the bipartition then 
the superoperator has trivial action. Therefore, the probability q(A) that 1Z{Ta) has a trivial action is given by the probability of 

drawing an edge completely inside A or completely inside B = V\A, that is, q(A) — K 2 ' , y N ~. 2 — -. Otherwise, TZ(Ta) has a 

(2) 

non trivial action: 

KW(T A ) - q{A)T A + ^-{N B J2 T A \i + N A J2 T Auj] (26) 

' ' iGA jeB 

Notice that the superoperator TZ^ acts on the space of shift operators T A , That is, for any subset A C V, 1Z( 2 '(Ta) is an 
operator in Hilbert space 'H.® 2 . Therefore, the matrix elements of V} 2 ^ in the space of the shift operators are given by 

KbU = (Tb\kW \T a ) (27) 
Notice now that the k— th iteration of the superoperator corresponds to matrix multiplication, that is, 

(n^) BA = (T B \(n^) k \T A ) (28) 



In order to find a convenient expression for 7Z for this model, we show how to map the superoperator 1Z(Ta) in Eq.(26 1 acting 
on the swap operators to a spin operator of the 2 w -dim Hilbert space (we are assuming now we are dealing with qubits). First 
of all, consider the mapping between subsets of V to a pure vector which is a member of the computational basis in the 2 w -dim 
Hilbert space of the qubits: 

K : A C V -> ® iev \ XA (i)) = \ip A ) 

wi*xA(i)={;;^ (^) 

Here, Xa(^) i s the indicator function for the node i £ V. By considering all the subsets X in the power set V(V) it is obvious 
that the states \ipx) , X C V form a complete orthonormal basis for the 2 N dimensional Hilbert space of N spin 1/2 particles. 
We will call this space the abstract qubit space and operators in this space are denoted with a hat on top. On this space let us 
introduce the total spin operators S a = g X^gy G t wnere erf are the Pauli matrices acting at the i'th spin. Note that then: 

S Z \^ A ) = \{N A -N B )\^ A ) (30) 
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FIG. 1. [i] (left panel) Purity as given by Eq.i 38 i with fixed total size N 
2 (purple squares), N A 



10 and different subsystem sizes N A = l(blue circles), Ljv = 
4(green upright triangles), Na ~ 5 (dark blue inverted triangles) [ii] (right panel) loga- 



3 (yellow rhombi), Na 

rithm of the difference between purity at the fc'th step and the asymptotic purity with N = 10 and same color codes for Na as in [i 



which allows us to promote the numbers N A , N B in Eq.(26 1 to operators Na —¥ Na = (S z + N/2),N B -> N B = {N/2 - S z ), 



N A N B 



1 



in the sense that N A — {ip a\N a\i}j a) and similarly for N B . With these definitions, we have qa = 1 — n(n-i) 
|-gy [(A^/2) 2 — (S z ) 2 ]. Similarly we now define the total raising S + = S x + iS y and lowering operators S~ = S x — iS y . Then 
we obtain 



S \^a) = ^2 \ipAuj) 



(31) 
(32) 



ieA 



Now we can define the (hatted) operator TZ^ acting on the abstract qubit space as: 



ft® := [(iV/2) 2 -(^) 2 ]) 

+ ^[S+(N/2 ~S Z ) + S-(N/2 + S z )} 
\hi\ 



(33) 



is a non-Hermitian operator in the abstract qubit space spanned by the vectors \i/j a )- In this space, its matrix elements read 
(RJW) b ,a = (V'bI T^ 2 ' IV'a)- We now show that Iterations of the RQC protocol in the original space correspond to iterations of 
the operator p3) l in the abstract space. To this aim, an important observation is that in the original space the trace of the powers 
of the superoperator 1Z is taken w.r.t. w® 2 which belongs to the totally symmetric subspace of "H® 2 . We can finally relate the 
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average purity of the subsystem A to a sum over matrix elements of the matrix (TZ^) k : 

* k (A)=Tr[u® 2 n ( k 2) (T® 2 )} 

E {K { ?)b,aT b ) 



P 



a=2. 



,®2 



BCV 
BCV 

= e (M& w ) k m 

BCV 
BCV 



and thus from Eq.( 34 1 one can see that in the abstract qubit space 

p^ k = m + {i\r N {n^) k \^ A ) 



(34) 



(35) 



because J2bcv IV'b) = + ll))®^ = IV'symm)- Notice that the state \ip s ymm) is invariant under the action of the projector 
Ilg = J2aes N a onto me s y mme tnc subspace of the abstract qubit space, where a € 5 at are the elements of the symmetric 



group on TV labels. With this in mind and the fact that Total spin operators commute with II5 we can write Eq.( 35 1 as: 

,„„, m |n s (^( 2 )) fe n s |^) 



p° 



(A 



(36) 



Note that II5 \iPa) € Hj=j max i.e. II5 is the projector onto the highest total spin subspace. In our case then it projects 
onto the J = N/2 subspace. Since [{TZ^) k , Us] = we focus on just the %j—j max subspace. However II5 \ip A ) is not 
normalized. Normalization involves some algebra to relate states in the tensor product basis for the abstract qubit space to the 
total spin basis: suppose then that the z-component of spin for some state in the tensor product basis of L spins is m. Then 



Ndown) = m 



N„ 



(N/2 + to) thus the number of distinct states in the tensor product 



basis with m as their z-component is ) = ( N /2+ m )- The state in the total spin basis in the %j-j max subspace that has the 
same z-component is a symmetric combination of these distinct states with appropriate normalization. It turns out that: 



n 5 |m) 



N \ 
N/2+mS 



\N/2,m) 



With this in mind we see that fig \ip/ 



into Eq.( 36 1 we obtain 



\N/2, \{2N a - N)). Inserting this and the identity in the H 



j=j„ 



(37) 



subspace 



N/2 

E 



\J (n a ) m=-N/2 

x (N/2,m\ (K {2) ) k 



U s \L/2,m) 



N/2,-(2N A -N) 



C(N,N A ) 



N/2 

E 

m=-N/2 



N 

m + N/2 



(7e (2) ) fe , 

y'^- Kn,\(2N A -N) 



(38) 



where C{N,N A ) = 
(^,m\(K^) k \"\(2N A 



N A \(N-N A )\ 
N\ 



and the matrix element of the k iterated superoperator are (TZ^) k 



AO). 



Note that Eq.(|38| expresses purity dynamics of a d N -dimensional system in terms of dynamics in an exponentially smaller 
(N + 1) -dimensional space. From eq.(33 1 one can see that for N A = 0, TZ^ \ip A ) = 1 as well as for N A = N, ti^ \^ A ) = l, 
all other eigenvalues of the operator being less than 1 . If one assumes that these are the only two fixed points of the operator then 
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log[A] 




log[LJ 



FIG. 2. [i] (Top panel) Scaling behavior of gap(A) for the superoperator of the Fully connected RQC. The best fit line(red) gives us the 
equation: log[A] = —.969395 log[L] + .0246366 [ii](Bottom panel) Scaling behavior of the logarithm of the product of operator norms for 
the similarity transform matrix M and it's Inverse. The best fit line(red) gives log[|| M\ \oa\ | A^" 1 1 1 oo] = 1.17508 + 0.318L 



asymptotically the non-zero eigenspace is spanned by the symmetric and antisymmetric combinations 
with only the symmetric combination contributing to the purity which reaches the value 



V2 ' 



\4>o)-\4>n) 

V2 



pc 



d 2N-N A + d N+N A 



d N (d 



N 



1) 



(39) 



The typical behavior of purity is shown in fig.fTJ. An obviously important question is : How fast in k does the protocol take the 
purity to within e > of the asymptotic value ? The answer to this question can be related to the gap of the matrix TZ^ in the 



subspace. From Eq.( 33 1 it is easy to see that the matrix representation R has elements: 



N 



K 



(2) 



N 



A ft \^ Nd U N 



+ (y-ff) 



N N 



In n 



9+1 



(40) 



which implies that the above matrix R is tri-diagonal and satisfies the condition Rk.k+i ■ Rk+i,k > 0. Such a matrix is 
similar to a Hermitian matrix. One can then diagonalize the Hermitian matrix and asymptotics can be calculated by the power 
method of eigenvalues for the derived Hermitian Matrix. We then have that SRS^ 1 = H 3 U s.t. U^HU = D => 

IP SRS~ 1 U = D where D,U, S are diagonal, Unitary and Invertible matrices respectively. Using the spectral resolution of D 
one can writcUSRS^W = J2i A * => Rh = S^t/ 1 E t >% (M us Then with M = c/t s we obtain from 
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Eq.d36j> for the case \A\ = N/2: 



N+l 

P^ k = C(N, N A ) £ \ k a <V + | M- 1 \X a ) (X a \ M 

a=l 
2 N/2 W +! 



' N * 
JV/2/ 



53 A « a « 



(41) 



ct=l 



where |?/; + ) is the normalized |f/; sl)mm ) and J?/; ) = |J = N/2,m = 0} are unit normalized vectors and a Q = | M \X a ) (X a \ M 
Assuming that the eigenvalues of R are so arranged that Ai = Aa = 1 and the rest are arranged in a non-increasing order we 



find: 



L+l 



\Pr 2 -P£T 2 \=C(N,N A )\Y,^ a \ 



a=3 



L+l 



< C(N,N A ) ^ |A Q | fc | <V + | M \X a ) (X a \ M- 1 \$°) 



a=3 



L + l 



< 



C(N,N A )\X 3 \ k 53 | (tp + \ M \X a ) || (A Q | M" 1 



Q = 3 



<C(iV,iV A )|A3| 



L+l 



53l(^+|M|A Q ) 

\ a=3 



L+l 



53l(A Q |Af^|^ ) 

\ a=3 



<C{N,N A )X^\\M\\ 00 \\M-% 



(42) 



We require the R.H.S of the above expression to be less than e. Therefore, taking the logarithm of the inequality C(N, N A )X k \\M\ | c 
e we obtain (using A3 = 1 — A) : 



k > 



Since log — ^> requiring 



k > 



logC(A^)+log||M|| 00 ||M- 1 |U+-logl/e 
1o St^a 

logC(^,^ A )+-log||M|| 00 ||Af- 1 || 00 +-logl/e 
A 



(43) 



(44) 



makes sure that inequality Eq.(43 1 is also fulfilled. Indeed the bound in Eq.([44|> is a rather weak lower bound as for all graph 
sizes N that we studied numerically the asymptotic values for any size of the subsystem N A < N were reached much before 
the above bound. 

The main difficulty in estimating k min (e, N, N A ), i.e. the minimum number k of iterations required to reach within e accuracy 
of the asymptotic value resides in the calculation of the operator norms of the similarity transform matrices M, its inverse M^ 1 
and the gap A. This amounts to diagonalizing the non-Hermitian matrix R in the maximal total spin subspace. However this is 
an exponentially reduced problem of diagonalization in a (N + 1)-Dim space compared to a d N dimensional one. Numerical 
study shows that the gap A of R to have an algebraic dependence on N as shown in Fig.Q. From the least-squares best fit line 
we can evaluate 



A = e- 025 /iV- 97 



1.025 

N 



(45) 



Moreover, from the lower panel of the same figure we find a linear dependence of the logarithm of the product of operator norms 
of the matrices M, M^ 1 i.e. 



log 



\M~ 



1.17508 + 0.3187V 



(46) 



Thus logC(£) ~l^oo ai logL, logfllMHooHAf 1 \\ 00 ] a 2 L + ft, ^ ^ a 3 L where ai,a 2 ,a a , /3 = 0(1) we finally obtain 
the scaling 



k > [an log N + a 2 N + /3)a 3 N + log (l/e)a 3 iV = 0[N 2 ) 



(47) 
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FIG. 3. A bipartite (^4, B) spin chain of length L = La + Lb with nearest-neighbor qubits interacting via 2— qubit gates (ellypses). The edge 
e is the one that straddles the two partitions. The gates are numbered by the subscript Xi where x = A, B denotes the two halves of the chain 
and i the distance from the boundary of the two partitions. 



V. CONTIGUOUS EDGE MODEL 



This model intends to mimic evolution of a multi-partite system under a local (time-dependent) Hamiltonian. Here, all the 
edges of graph Y are acted on by the RQC with 2-local unitaries in some particular order 'a' - denoting an ordered sequence 
of edges. A RQC of depth O(N) picks this ordering of nodes such that finally all nodes in the graph are acted on by nearest 
neighbour unitaries. One pass through such a circuit is called a cycle. This procedure is then iterated through n c cycles. The 
composite unitary of a cycle is the a— ordered product of unitaries with support on nearest neighbours according to the graph Y. 
The ensemble is constructed by considering all the possible orderings a. 

While our general formulation of this model extends to graphs with any geometry in any number of dimensions we present a 
detailed analysis of 1-D graphs and the 2-D square lattice. 



A. Contiguous Edge Model on the Linear Chain 



Let us start with the 1-D graph of Fig.Q, with N = L, and introduce a bipartition into subsystems A and B of lengths 
La, Lb- The sites in A on the left of the boundary are labeled by %a — I a, ■••> La increasing towards the left, while the sites in 
B are labeled by ig — Is, Lb increasing towards the right. The edges on the chain are labeled by a, = ((i + 1)a, Ia)> an d 
the edges in B are labeled by bi — (is, (i + 1)b- The RQC chooses an ordering among all the edges in the chain, and then on 
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each edge acts with a Haar distributed 2-qudit unitary operator in the given order. The ensemble therefore contains all possible 
permutations of the list of edges. For n c = 0(1), the average purity depends strongly on the order in which the RQC chooses 
the edges, while when n c exceeds the subsystem size, ordering does not really count, as we shall see in the following. In any 
case, we can consider the orderings that give the two extreme situations, that is, the minimum and maximum decrease of purity, 
i.e. maximum entangling power which we call the best sequence and the one that corresponds to the minimal decrease of purity 
or minimum entangling power is termed the worst sequence. Let U a denote the ordered product of 2-qudit unitaries over all the 
edges in E with the order given by the permutation a i.e. U a = U a { ei ) ■ ■ ■ ?7 CT ( e|E| ). The ensemble and its measure are then given 
by 

^(^ = {^1*)}^ 

dn(U) = 6(U - U a ) [] dn Haar {U e ) (48) 

We can now see why sequences of unitaries corresponding to different permutations yield different purities. We start showing 
which sequence produces the greatest decrease of purity, that we dub best case. Consider the sequence Ubest = U^UaUb where 
U e is the unitary straddling the edge i.e. acting on 1a and 1b and Ua = U ai U a2 ...U aLA _ 1 is the internal structure of Ua 
where U ai means a 2-qudit unitary with support on the qudits the nearest of which is at a distance of 1 lattice spacing from the 
boundary. Note that the internal structure of the unitary is non-deformable since [U ai , U ai+1 ] ^ as they share a node. With 
the same convention the internal structure of Ub — Ub 1 Ub 2 ...Uf )Lg _ 1 . Physically this corresponds to the RQC choosing and 
applying all possible unitaries on the outermost nodes in A and B followed by two nodes 1 lattice spacing closer to the boundary 
and so on till the nodes 1^, 2a and Is, 2 B are acted on. Finally there is a boundary interaction through U e . In terms of the 
probability distribution according to which the edges are picked by the RQC the best sequence is a cycle of length(circuit depth) 
k = La + Lb — 1 where the (La — 1) levels of the circuit choose edges in A as follows: The first set of nodes X\ is chosen as 



V {1) (X 1 ) 



liff|Xi|=2 A D(X u l A ) = L A -2 
= otherwise 



(49) 

while the next L A ~ 2 levels choose edges depending on the choice of nodes in the previous step 

1 iff \Xi_i HXi\ = 1 A 

D{X u l A ) = LA-{i + l) 
= otherwise 

(50) 



V^(X i \X i _ 1 ) = I 



where i = 2,3, (La — 1) and the distance between two sets D(Xi,Xj) = min x (\xi — xj\) is the minimum difference 
between any two elements belonging to the two different sets. In our case this difference is the difference of the position labels 
(1, 2, La). Similarly of the next (Lb — 1) levels of the circuit the first set of nodes Xl a in B is picked as follows: 

v {La \x La \x La ^) 

= fl iff \X La \ =2 A D(X La ,1 b ) = L b -2 
1 = otherwise 

(51) 

while the next (Lb — 2) sets are chosen as 

f 1 iff \Xj n Xj^l = 1 A D(X h 1b) = L b - (j + 1) 
1 = otherwise 

(52) 

where j = 2, 3, (Lb — 1). Finally the RQC chooses the boundary edge as follows: 

V {La+Lb-1) {Xla+Lb _ AXla+Lb _ 2) 

= fl iff \X La+Lb ^\ =2 A \X LA+LB - 1 r\A,B\ = 1 
1 = otherwise 



(53) 
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One can similarly devise probability distributions for the RQC that generates any desired ordering of the edges. However, for 
now, with the measure defined in Eq.(|48]l and iterating the algebra of Eq.([T4|, we find that the purity is given by nested integrals: 



P«=^TY[Lo® 2 Jd[U B ]d[U A ]dU e {ulu\uiT 2 T A (U e U A U B )® 2 \ 
= Tr[c^ 2 / d[U B ]((U B r 2 I d[U A ]((U{)® 2 



( / dU e (Ul)® 2 T A Uf 2 )Uf)U® 2 )] 



= Tr[oj® 2 (N 2 T A + JVjT A _i + iV d 4 T A _ 2 .. + N^ a T a _ La+2 
+ N^T a _ La + N 2 T A + NjT A+1 + N*T A+2 ... 

+ n^t a+Lb _ 2 + n^t a+Lb )} 

= (N 2 + N$ + Nj + .... + + N^)+ 

(N 2 + Nl + Nj + .... + N%» + N^) 

_ N 2 (l-N^) {La _ 1} 
(l-N d ) + d 



N^l-N^- 1 ) , (i, B _i) 



(1 - N d ) 

(54) 

The notation T A+r means the swap acting on X = A U 1b, %b, tb and similarly X = A — r = A\1 A , 2 A , r A . For 



L A , L B reasonably large we find from eq.(54 1 that at the conclusion of the first cycle i.e. iteration 1, 



while for n c <C L A , Lb 



OjV 2 

»c=l~ il _ Nd y (55) 



P^n^lT^L) , (56) 



(1 - N d ) 



and again we recall that this is for the best case sequence. Also note that from Eq.( 14 1 we see that for each non-trivial action of 
the averaging procedure over the unitaries we get a decrease of purity by the same amount,(l — 2N d ), thus a sequence which 
maximizes the number of non-trivial actions will result in the maximum decrease of purity. 

Now we want to compare the above case with the sequence that produces the least decrease in purity (or the least 2-Renyi 
entropy), or the worst case. We see that by keeping the same internal structure of U A , Ub and comparing the purities corre- 
sponding to the four possible cases: (l)U = UJJ a Ubi (2)C/ = UbU s U a , (3)C/ = UjJJJJb, (4)C7 = U B U A U e we can easily 
see that for n c = 1 the sequence (4) is indeed the worst case scenario. Numerically we find that sequence (4) performs the 
worst also for generic large n c . As for n c > 1, the decrease of purity also depends on what is the ordering of the unitaries in 
the products U A and Ub- We find numerically that the least decreasing sequence of unitaries (worst case) is given by choosing 
Uworst = U B U A U e with U A = U aLA _ 1 ....U a2 U ai and U B = U bLg _ 1 ....U b2 U bl . Let us quantify the purity using the worst 
sequence U wors t as function of n c . For n c — 1, we have 



P a=2 n c =i = T^fa (NctT A _i + N d T A -i)] = 2N d (57) 

We can also obtain the exact expression for the purity after any number of iterations n c < L A where the size of the environment 
L B > L A . We find: 

pa=2 

n c 

Tit — Th c X / \ 

= J2 TrK^^|^^ +m )(T A _ (m+1) +r A+(m+1) )] (58) 

m=0 

= 2s n <{l - ((1 - N d )N d ) n ° {^ nc ~ l \ 2 F 1 (1, 2n c - 1 + n c ; N)} 

2N n a 

d -{l-f(n c )} (59) 



(i - N d y 
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where C(n c , m) — 2( na+ ™ 1 ) , s = n~^7J ' anc ^ 2-Fi(l, 2n c ; 1 + n c ; 2V) is the Gauss Hypergeometric function. Therefore, our 
formalism allows us to obtain typicality of the purity for arbitrary depth of the RQC. When n c — O(l), the subsystem is very far 
from being the Haar case, and the RQC is not a t— design for any t, indeed, the system features an area law for the entanglement. 
As n c increases, entanglement propagates in the bulk at distance ~ n c from the boundary. For large values of n c , the expression 
simplifies because f(n c — > Large) —> 0, and we see that 



2/V™° 

pa=2 pa d / 6 q^ 

Comparing this expression to the one for the purity in the best case sequence Eq.d56k we see that, as long as n c is smaller than 
the system size, the best case sequence is better by a factor of i(2A^)" c , see Figl4 top panel). Nevertheless, the same figure 
shows that after system size is reached, the two cases converge to a similar value. In the following, we show that the asymptotic 
value is independent of the ordering. This means that there is a n c above which there is an onset for the independence of the 
ordering. The numerical results shown in Fig.Qtop panel) suggest that the onset happens at n c ~ La- 

We can give a justification of why the above given sequences are indeed the ones that decrease the purity the most (best) 
or the least (worst). At n c = 1, of course, all the unitaries acting after U e are not entangling at all. With n c growing, all the 
unitaries that are acting before U e allow some entanglement to be generated among the qubits. The sequences that go towards 
the boundary bring entanglement towards it, while the ones which start with the boundary and go outwards bring entanglement 
away. In this case, indeed the site 1a would get very entangled with the bulk in A (the same occurs to the B side), and by 
monogamy this does not allow to effectively transfer entanglement information across the partition. This scenario also shows 
how in the worst case the RQC would entangle nodes at the same length as the iteration number. 

At this point, we want to look for a result about the maximum purity of the n c — iterated ensemble £" c . The purity is a 
monotonically decreasing function because it is obtained by iterated application of CP maps, for all n c . As we pointed out 
above, the different results above hold in the region of validity for the scaling of n c . Of course the absolute minimum of the 
purity of the reduced state cannot be less than (^) La corresponding to the totally mixed state. To answer questions about the 
average distance of a state in £ n " asymptotically to the totally mixed state on the subsytem A we resort to the superoperator 
formulation a little later. As we noticed above, plot in Fig.(|4jTop panel) also shows that for large n c the ordering does not count. 
This means that our scheme does indeed mimic the Trotter scheme as far as the statistics of the reduced system is concerned. In 
other words, the average decrease of purity in this model approximates the average decrease of purity obtained by evolving with 
local time dependent Hamiltonians and using the Trotter scheme. The irrelevance of ordering for large n c can be understood 
mathematically in the superoperator formalism, what we now go on describing for this model. The chain superoperator is a 

a— ordered product of non commuting projections i.e. lZ c hain = T^a(\E\) 7£o-(i)> because of the non-commutavity of the local 

unitaries that make up the product. Infact the products in the superoperator is in reverse order of unitaries. For the sake of brevity 
without explicitly writing a = 2, using eq.|9]) and the first line of eq.(54i, one can see: 



U = U n U n -i....Ui =^> TZchain = 7?-i....72. n _i72.„ (61) 



where the subscripts on U denote supports for the same. Note that, for a fixed sequence, V} chain ^ Tt c hain and therefore lZ c hain 
is not a hermitian operator. Nevertheless, the averaged sum over all possible sequences, 7Z — X TZ a , is Hermitian. Since 
our definition of the best sequence is U — UJJaU b with a specific internal structure of unitaries within Ua and Ub and that for 



the worst sequence is U = U A U B U e Eq.(61 1 implies that IZbest = T^ worst . One can understand the action of the superoperator 
by studying it's action on the non-Orthnormal basis of the swap Tx, X C V operators. We regard the swap operators as kets 
in the + 1)— dimensional subspace S = span{|i)}, < i < \V\ of the Hilbert-schmidt space on (Wy) 82 . This subspace 
is the space of swap operators acting on all i qubits from one end of the chain where < i < \V\. The correspondence 
\i) = Ti = i-vu ® Ti then implies that the Hilbert-schmidt inner product ^ Sij. In this basis the matrix representations of 
these superoperators are real and so it turns out that the eigenvalues of lZt, es t and that of lZ worst are identical. This means that 
for sufficiently large iteration number n c which sequence we choose does not matter while the rate of approach to the asymptotic 
value of purity is dictated by just the gap (1 — A2) in either case. The difference in initial decays of the purity for the sequences 
lies in the fact that the eigenvectors correspoding to identical eigenvalues are different. 

Let us now explain the action of lZ c h a i n m our chosen basis. Consider the action of the superoperator TZ c h a in corresponding 
to the best sequence for the linear chain which has \A\ = La, \B\ = Lb when it takes as an argument some Tx , X C V which 
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FIG. 4. [i](Top panel) Purity as a function of n c for the best sequence of unitaries(circles) and for the worst (squares) for La — Lb = 8. 
The smooth line connecting the squares is the analytic expression for worst case P„ c . [ii] (Bottom panel) Convergence of the worst case 
sequence to the asymptotic formula value for different values of system/environment length. Blue circles(L = 5, La = 2) and Yellow rhombi 
(L = 6, La = 3) converge to > 99% of their asymptotic value whereas Purple squares(L = 8, La = 3) converges to > 97%. 
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is the swap on (Hx)® 2 subspace. Then the possibilities are : 

n chain (\o)) = i, n cham (\v)) = i 

n chain (\i)) = N d +P I* + 1 - P) + N d |0) , K * < La 

p=0 

l a -i 

K chain {\L A )) = J2 N 1 d +p \L A + l-p)+N^\0) 
P =i 

\V\-La-1 

+ Yl N 1 d +p \L A -l+p) + N^\2L A ) 
P =i 

\V\-i-l 

n chain (\i))= N l +p \ i - 1 +p) + N d Vl ~ i \ v ) 

p=0 

,L A <i<\V\ 



(62) 



In this basis the matrix representation R of H c hain cannot be guaranteed to be even Normal thus we may only attempt a Jordan 
decomposition of R but this is enough to understand the iterative behavior of the RQC. The purity of the generated ensemble 
after any number of iterations k is given by: 

P^n c = ]T (i\ R n ° \L A ) (63) 
\i)es 



From Eq.(61 1 we see that ||72- c hain|| < IleeB ll^ell < 1 which means that all eigenvalues A, are less than equal to 1 in 
modulus and hence asymptotically only the contribution from fixed points survive. Note that each of the projections lZ e has two 
eigenvectors 1 and TV (which are swaps on no nodes and all nodes respectively), with eigenvalue 1 . If one assumes then that 
the common eigenspace of TZ c h a i n is spanned by 1 and Ty we find that the asymptotic value of the purity is: 

P^ nc ^ 00 = {u®\K™ hain (T A )) 

V ' y/d L (d L + l) n y/d L (d L + iy ' 



/ ®2 t-Ty t-Ty 

^ ' y/d L (d L + l) n y/d L (d L + l)' ' 

1 ' 1+TV T A) 



^d L (d L + l)\/d L (d L + 1) 

d 2L-L A + d L+L A 

dW + l) (64) 



where the normalization for the symmetric and antisymmetric combinations of the basis vectors for the common eigenspace is 
obtained by setting ( ]L +^ V ' , = Tr[ ^^ ] = 1 =^ C = ^Jd L (d L + 1). We have verified this result numerically using 

small system lengths and the convergence with n c to the value given by eq.([64]) is shown by fig.Q. 

For the specific case of L A — Lb we found numerically that the sub-dominant eigenvalue A2 saturates at a value of (2Nd) 2 
Pig.(j5]» with increasing system size. Indeed numerical evidence indicates that the spectrum of TZ c hain at l east f° r this case is 
entirely positive < A : A 6 Spec(7?. c / la i„) < 1. All that happens upon increasing the system size is that the population of 
eigenvalues in the subregions of the domain < A < A2 increases proportionally, except for the largest eigenvalue of 1 whose 
number remains equal to 2, see Fig.|5]l. This corroborates our assumption that there exist only two fixed points. With this value 
of the asymptotic purity we can bound the average distance of states in our ensemble reduced to the subsystem, A, from the 

totally mixed state on it, u = ■ 
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FIG. 5. [i]Saturation of the sub-dominant eigenvalue for La — Lb with increasing La for a chain of qubits. Similar behavior is observed for 
general qudits.[ii]Distribution of eigenvalues for a chain of qubits with La = Lb- The orange bars represent population levels within bands of 
the domains for the eigenvalues with La — Lb = 200. Blue bars(superimposed on orange) are for La = Lb = 100 The gap in the spectrum 
can be clearly seen as the difference between the largest eigenvalue of 1 and the next =.64 
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Indeed we see that the trace distance, V(p,a) — \ ||p — cr| |i < \ \J 'rank(p — &)\\p — ex 1 1 2 while 



cr||p -oj| 2 = - oV Tr [P 2 ~ pcr-crp + cr 2 ] 





P{p) - P(a) 



(65) 



Thus for S B \ S A)[ we can get indistinguishably close to the totally mixed state of the subsystem. 

The analytical calculation of the variance for this model is difficult because the unitaries do not have disjoint supports. How- 
ever, we can resort to a Markov type inequality for positive valued random variable to assert that the ensemble has a small 
variance since the average purity itself for large number of iterations reaches a value exponentially small in the system size. 

To summarize, the ID model shows that for n c = 0(1) almost every time dependent local Hamiltonian will evolve a factorized 
state to a state obeying the area law. The fluctuations around this average are small in view of the Markov inequality. On the 
other hand, when the number of iterations (the "time" in the Trotter scheme) scales with the size of the subsystem n c > La, 
almost every evolution will produce a state with the volume law, and, asymptotically, maximally entangled. 



B. Contiguous Edge Model on the Square Lattice 

Of course the area law in ID is rather trivial. In order to find a more meaningful result, we need to study a two dimensional 
situation. We will see that in this case the fluctuations around the area law are even smaller. We consider a 2-D square lattice 
with a bipartition into A U B as in Fig|6] In this case the boundary between A and B is obviously a Id system. For sake 
of simplicity, let us start with a RQC that applies Haar distributed 2-body gates across nodes of the bipartition that are nearest 
neighbors. It is easy then to see that the unitaries across the boundary will have disjoint spaces of action and the average purity 
is a product of the single edge purities for the n c = 1 iterated ensemble. Explicitly , 

= Tr[^ 2 | dU ei (Ul f 2 T lal U® 2 J dU e2 {U\ 2 f 2 T la2 U® 2 

... J dU ei {Ul i )® 2 Ti al Uf 2 T A \ lal la2 ,,^i a i\ 
= Tr[uj m N d (li a iAbi + T la iTi b i).N d (li a2t ib2 + T la2 Ti b2 ) 

■ ■■■N d (\i a i :1 u + Ti a iTi b i)T A \i 

al.la2,..,la/J 
= Tl[u (g>2 N l d (li a i j i b i^ a 2 : ib2,....lal,lbl + 
1^161.1 a2,162, la/, lbl T" ■■ 

al,la2,..,la/J 

= (2N d ) 1 (66) 

where I is the number of closest nodes on either side of the boundary i.e. the length. This means that the 2-Renyi Entropy for 
this case can be lower bounded as S 2 > —\og(P a=2 ) — —I \og(2Nd). In words: the entropy is greater than a constant times 
the length of the boundary.In this case too it turns out that acting with the internal unitaries first we can get lower purities 
and the reason is the same as the monogamy arguments presented for the linear chain analysis. For eg. in Fig.(|6]l if we 
consider the very simple system of just the nodes l a i and \ a2 in A and l^i and \ b2 in B then, upon implementation of the 
sequence U = Ui b i.i b2 Ui a i.i a 2Ui a i,i b iUi a2t i b2 the average at iteration 1 is P a=2 k=i = (2Nd) 2 = .64 while for the sequence 
U = U lahlbl U la2tlb2 U lbhlb2 U lalila2 it is P a =h=i = 2Nj + 8Nj = .525. 

Let us now compute the variance. For the 2-D model we find that since the edges have disjoint supports the average of the 
square of purity is a product of the average of squares of purity for each edge and similarly the average purity itself is a product 
of the average purity for each edge: 

(P a=2 )l D = (P a=2 )l {P a - 2 )l 2 ....{P a = 2 ) 2 ei (67) 
PW =P^P^....P^ (68) 
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FIG. 6. A bipartition (A, B) of a square lattice. 

and the variance is thus 

Var[P^ 2 } =I l - (2N d ) 21 

( 2(2rf 4 + 9ri 2 + l) V / 2d V 1 

\{<P + Z){<P + 2)(<P + 1)) \(d 2 + l)J " ' 

which, if we have qubits d = 2 gives an exponentially decaying variance of 

Var[P^ 2 ] = -.64' + .657* (70) 

This is a strong result. The average 2-Renyi entropy does follow the area law. We are averaging over all the possible states 
obtained by local unitary transformations starting from a completely factorized state. This ensemble contains all the ground states 
of local Hamiltonian (without topological order). Moreover, we have shown that deviations from the area law are exponentially 
suppressed. Proving the area law in 2D for the ground states of local Hamiltonians is one of the most sought after results in 
quantum many body theory and our result represents a progress also in this direction. 
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VI. CONCLUSION AND DISCUSSION 



The question whether a physical model can generate an ensemble of states that reproduce moments over the Haar measure, 
that is, of the set of all the states in the Hilbert space, is an important question in quantum information theory, simulability of 
quantum many-body systems, and, recently the foundations of quantum statistical mechanics. In this paper we present some 
physical models with the aim of studying the statistical moments of the reduced subsystem A of an initially bipartite quantum 
many-body system A U B, and with locality constraints that should contain the essence of the local dynamics induced by a local 
Hamiltonian. In this formulation, the relevant size La at play is that of the subsystem A. 

We present a superoperator theory that allows us to make statistical claims beyond the Haar measure. In particular, we are 
able to discuss both the typicality of area law (which is quite far from the Haar measure) and the volume law (which, indeed is 
what would hold for generic states in the Hilbert space). In simple words, we can show that states evolved with a local evolution 
for times of 0(1) with respect to La, exhibit typicality of the area law for the entanglement, while for large times scaling with 
the subsystem size they exhibit a typical volume law. We want to stress the fact that typicality means that in the ensemble 
considered, the average is such and that the variance is vanishingly small in the large system limit. In other words, almost all 
members of the ensembles constructed show the above mentioned entanglement properties. In the asymptotic case, we recover 
results similar to those of t— designs, but again, where the scaling is that of the subsystem and not of the whole system. We have 
also discussed the implications of these findings for the foundations of statistical mechanics. 

Our results are based on the algebra of the swap operators on subsets of the total nodes in the graph. This algebra shows 
how the underlying graph-theoretic structure of the system determines the propagation of entanglement within the subsystem. 
In particular, the mixing time for the subsystem is strongly dependent on the connectivity of the graph and other graph-theoretic 
notions like the average distance between two nodes or the Hausdorff dimension. Recently, it has been advanced that black holes 
are fast scramblers, that is, systems in which the mixing time is maximal [43|,[44], and several models for the fast scrambling 
process have been put forward [43 ]-[47 ]. We believe that our techniques can be useful also for this line of research. 

As we discussed in the introduction, as long as topological order is absent, ground states of local gapped Hamiltonians can 
be obtained by a circuit of fixed depth from a completely factorizable state. It would be interesting to study the statistics of 
the entanglement in ensembles where the fiducial state is topologically ordered. For instance, we would like to know if in such 
ensembles there is a non vanishing topological entropy |48 1-|50| on average, and what are the fluctuations. This technique may 
prove useful to study the problem of the stability of topological phases. 
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For the symmetric group S a of order a, the elements are the permutation operators a on some space H where H can be 
broken up into a direct sum of Irreducible subspaces denoted by A each with dimension d\ and multiplicity m\ i.e. formally 
T-L = S>\£iRR(s a )G mx <E> C dx . For any operator a on the same space we have that: 
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Appendix A: Calculation of product of traces usiing Schur's Orthogonality Theorem 




(Al) 




(A2) 
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Using the above 2 equations we get that: 

l^Tr^Trta) 

a 

= cj E E to > to ^' x E E^fM^ 

ffeS„ \,X'eIRR(S a ) i,j=l 1=1 



E m ^'E(s E ^WV 

A, A' ij'.Z <t6S„ 



2 d * 



E^E4 A) = E^(^) (A3) 

A i— 1 A 



where in the second line above the Kronecker deltas appear due to Schur's great Orthogonality theorem [40|. 



Appendix B: Calculation of Variance 



The trace w.r.t. II + in Eq.(17i can be given a closed form expression again by using Schur's Orthogonality relation for 



representative functions on the irreps of a group see for eg.|40|). The result is 

Tr[n + (lf 2 > ® if 4 >) = J2 ^ T r[(7f 2) ® T^)U X }. (Bl) 

A:IRRsofS 4 A 

Notice that now the sum is over just the irreps of S4. The symmetric group of order 4, S4 has 4! = 24 elements which are all 
permutation operators a on 4 labels. In our case the labels refer to copies of the i'th or j'th space that make up (Hi ® "Hj)® 4 
. Thus a = (Ti ® <jj is a valid decomposition of the permutation operators. The projector onto the totally symmetric part of 
(Hi ® 'Hj)® 4 takes the form IT + = ^ ^2 aeSi cr, ® crj and hence 

Tr tJ [(lf 2) ®lf 4) )n + ] 



1 £ ^[(if a > ® if 4 ) ® ® a,)] 
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= ^ E ^tm (12) ®^.«7 i )].Tr J -[(l J , C r J -)] (B2) 

Note that in cycle notation the operator T^ 12 ' ® y^ 34 ) _ (12)^(34)^. Thus it's product with any other Oi still gives us a 
permutation operator on the % space. One can then just count the number of cycles in the product obtained where each cycle 
contributes a multiplicative factor of the dimension d of the space. For e.g Tr,-(<Tj = (12)^(34)^) = d 2 , Tri((7.j = lj) = d 4 since 
there are two cycles in the permutation (12)j(34), and 4 in Similarly the traces over the jth spaces can also be obtained. 
Finally the calculation above yields Tr id [(T$ 12) ® jf 4) )II + ] = d2 ( 2d4 + 9rf2 + 1 ) 
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